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Radio emission from cosmic ray air showers has the potential to become an additional, 
cost-effective observing technique for cosmic ray research, being largely complementary 
to the well-established particle detector and air fluorescence techniques. We present 
Monte Carlo simulatios of radio emission from extensive air showers in the scheme of 
coherent geosynchrotron radiation from electron-positron pairs gyrating in the earth's 
magnetic field. Preliminary results of our simulations are the predicted frequency, pri- 
mary particle energy, shower zenith angle, shower azimuth angle and polarization depen- 
dence of the radio emission. These properties can be directly related to data measured 
by LOPES and other experiments. 
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1. Introduction 

Accompanying the experimental efforts of the LOPES projeclpl, we are modehng 
the radio emission from cosmic ray air showers in the scheme of coherent geosyn- 
chrotron radiatiorP, thereby building a theoretical foundation for the interpretation 
of experimental data. As a follow-up to our analytic analysis of the emission proper- 
ties and the underlying coherence effectJ^, we have carried out sophisticated Monte 
Carlo simulations based on analytic parametrizations of the air shower proper- 
ties to achieve a high-precision modeling of the emission including its polarization 
characteristicJ^. Here, we present some of our preliminary results. 

2. Monte Carlo code 

Our Monte Carlo code calculates the geosynchrotron radiation from the air shower 
cascade's electron-positron pairs gyrating in the earth's magnetic field. All calcu- 
lations are done in the time-domain, without the application of any far-field ap- 
proximations. The full polarization information of the emission is preserved. This 
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approach is robust, fast and flexible and has been thorougUy teste(3. The underly- 
ing air shower model is based on analytic parametrizationg^ and takes into account 
longitudinal and lateral particle distributions, particle track length and energy dis- 
tributions, the air shower and magnetic field geometry and the air shower evolution 
as a whole. This allows a direct comparison of the analytic and Monte Carlo results. 
At a later time, the air shower model will be substituted by a direct interfacing to 
the air shower simulation code CORSIKA. Additional emission mechanisms (e.g., 
Askaryan-type Cerenkov-radiation) can be included at a later time. 

3. Simulation Results 

In this section, we present a number of preliminary results of our simulations. 

3.1. Frequency Spectra 

The spectral dependence of the emission (not shown here) clearly favors low frequen- 
cies for the measurement of radio emission from EAS. At a frequency of 10 MHz, 
the emission stays coherent up to high distances (~ 1,000 m) from the shower 
centre. Additionally, the signal strength is higher than at higher frequencies. At 
a frequency of 55 MHz, coherence holds up to ^ 300 m for a vertical air shower. 
Observing EAS at larger distances, i.e. in the incoherent regime, is still likely to be 
possible, but Monte Carlo simulations with a more sophisticated air shower model 
as will be achieved through interfacing with CORSIKA are needed. 

3.2. Primary particle energy 

In the coherent frequency regime, e.g. at 10 MHz, the electric field strength scales 
almost linearly with the primary particle energy (not shown here). This is the 
expected behavior for coherent emission. If an increase in depth of shower maximum 
as a function of primary particle energy is taken into account, the energy-scaling 
becomes dependent on the distance from the shower maximum. While it stays a 
power-law at all times, the power-law index drops from ^ 1.1 in the centre region 
to ~ 0.8 at 500 m distance. 

3.3. Polarization characteristics 

Figure^) shows the total field strength contours for a 45° zenith angle air shower. 
The emission pattern is elliptic due to projection effects arising from the inclination 
of the shower axis. Overlaid are indicators denoting the ratio between east-west and 
north-south linear polarization. As the total field strength pattern is simply aligned 
with the air shower axis and does not reveal any significant dependence on the 
geomagnetic field direction, it cannot be used for a verification of the geomagnetic 
emission mechanism. However, the emission is highly polarized in the direction 
perpendicular to the shower and geomagnetic field axes. Experiments measuring 
linear polarization can thus directly test the hypothesis of geomagnetic emission. 
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Fig. 1. a) Total 10 MHz field strength contours for a 10^^ eV air shower with 45° zenith angle 
coming from the east, overlaid with polarization indicators denoting the ratio between east-west 
and north-south linear polarization, b) Zenith angle dependence of the 10 MHz emission. 

3.4. Zenith angle 

As can be seen in Fig. ^p), the radial dependence of the emission flattens sig- 
nificantly with increasing zenith angle. Apart from simple projection effects, an 
intrinsic flattening occurs as a consequence of the growing distance of the shower 
maximum for increasing zenith angles. The slight deviations from the trend visi- 
ble for the 15° case (i.e., at only 5° relative to the 70° inclined geomagnetic field) 
demonstrates the weak influence of the geomagnetic field on the field strength. 

4. Conclusions 

We have developed a powerful Monte Carlo code for the simulation of radio emis- 
sion from cosmic ray air showers. Our predicted dependence of the radio emission 
on important air shower parameters provides a necessary prerequisite for cosmic 
ray research with radio techniques. In particular, the predicted polarization char- 
acteristics allow a direct verification of the geomagnetic origin of the radiation with 
LOPES. A detailed analysis of our results is in preparation. 
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